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Abstract 



The quantum-catastrophe (QC) benchmark Hamiltonians of paper I (M. Znojil, J. Phys. 
A: Math. Theor. 45 (2012) 444036) are reconsidered, with the infinitesimal QC distance A 
replaced by the total time r = vl — A G (0, 1) of the fall into the singularity. Our amended 
model becomes unique, describing the complete QC history as initiated by a Hermitian 
and diagonalized A^— level oscillator Hamiltonian at r = 0. In the limit r — ?■ 1 the system 
finally collapses into the completely (i.e., A^— times) degenerate QC state. The closed and 
compact Hilbert-space metrics Q^^\t) are then calculated and displayed up to A^ = 7. The 
phenomenon of the QC collapse is finally attributed to the manifest time-dependence of 
the Hilbert space and, in particular, to the emergence and to the growth of its anisotropy. 
A quantitative measure of such a time-dependent anisotropy is found in the spread of the 
A^— plet of the eigenvalues (r) of the metric. Unexpectedly, the model appears exactly 
solvable - at any multiplicity A^, the A^— plet of eigenvalues (r) is obtained in closed 
form. 
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1 Introduction 



In paper I [T] we analyzed bound-state energies Er'{\), n = 0,1,...,N - 1 generated 
by certain toy- model Hamiltonians H^'^\X). A characteristic property of our model was 
that an infinitesimal decrease of the real parameter A below critical Aq implied an abrupt 
complexification of the whole spectrum. Such a multifurcation quantum catastrophe (QC) 
was interpreted as an A^— tuple analogue of the Thom's bifurcation catastrophe which often 
occurs in classical dynamical systems and which is usually called "cusp" [2]. 

The transition from classical to quantum catastrophes proved facilitated, first of all, by 
the recent new developments in quantum theory. A review of these developments may be 
found, e.g., in Ref. [2] (where the underlying innovative representations of the operators of 
observables were introduced and called quasi-Hermitian) or in Ref. [1] (where the author, 
very successfully, recommended the use of special, PT— symmetric and PCT— symmetric 
quasi-Hermitian Hamiltonians, with V representing parity, with C being a charge, and with 
symbol T denoting the usual time reversal) or in Refs. [5] (where one can find an immediate 
connection of the new physics with the traditional mathematical theory of operators which 
are self-adjoint in Krein space). 

In Appendix A of paper I the quasi- Hermiticity property of the Hamiltonians H^^\\) 
was re-christened into "crypto- Hermiticity" , replacing its widespread but rather misleading 
nickname of "non- Hermiticity" . A compact account of the terminology may be also found 
in Ref. [B] where we recommended to call the whole formalism the "three-Hilbert-space 
(THS) quantum mechanics". In our present paper we intend to make the next step in 
connecting such a THS formalism to the currently open problem of the definition and 
systematic description of catastrophic scenarios in quantum theory. 

Several topics will be addressed. Firstly, in section [2] we shall briefiy summarize the 
preceding developments and recall some properties of our benchmark bounded-operator 
Hamiltonians. Next, in section[3]we shall discuss the necessary specification of the suitable 
physical Hilbert space, paying particular attention to the ambiguity of the definition of the 
mathematically correct inner product in this space. Our present innovative, QC-adapted 
recipe of the removal of this ambiguity will be thoroughly analyzed in section |4] where we 
shall present arguments in favor of the choice of parameters which would, in some sense, 
minimize the anisotropy of the physical Hilbert-space metric 0. We shall also replace there 
the time t of paper I (which starts at the QC singularity) by the more natural time r which 
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runs in opposite direction and which is assumed to have its initial zero long before the fall 
of the system into the final and singular QC degeneracy. 

In section we then turn attention to the systems with the smallest level- multiplicities 
< 4. Using the brute-force linear algebra methods we shall construct the fully explicit 
forms of the matrices of the metrics which we declared, in the preceding text, optimal 
and unique. On this background, the qualitative features of the QC process will be shown 
related to the properties of our manifestly time-dependent physical Hilbert spaces. 

In section [6] the transparency and compact form of the previous results and formulae 
will enable us to employ extrapolation techniques. Their use will extend the validity of 
many of our previous considerations to arbitrary A^. The whole picture of the A^— tuple 
QC level-degeneracy scenario will be finally summarized in section [71 

2 Benchmark models 
2.1 Complexified potentials 

In the current literature the applicability of the THS formalism is most often illustrated 
by the replacement of the linear harmonic oscillator Hamiltonian H^^^'^^ = + by 
the Bender's and Boettcher's [7] set of anharmonic models H^^^\5) = + x^^^ . They 
are distinguished by the real exponent 5 > 0. In a historical perspective their study 
appeared well motivated by the needs of quantum field theory (cf., e.g., [8]) and of various 
perturbation expansion methods [9| [TO]. 

In contrast to the exceptional selfadjoint operator if^^^'^) = H^^^\0) which lives in 
the most common representation L^(M) of the physical Hilbert space of states, the correct 
physical interpretation of the apparently non-Hermitian Hamiltonian H^^^^ (S) requires an 
ad hoc specification of the correct physical Hilbert space of states 7i^^^ ^ L^(M) In 
other words, the correct probabilistic quantum predictions are only obtained when one 
redefines the inner product in the vector space of the square-integrable wave functions in 
a way explained in Refs. |3l [6] , 

(V'll^a) ^ (Vi|^2)(s):=(^i|e(^)|^2). (1) 

Needless to add that the metric 0'-'^) ^ I must satisfy severe mathematical restrictions as 
listed in Refs. [21 [TT]. Only then, the standard physical Hilbert space H^^^ becomes well 
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defined and the Hamiltonian itself may be declared self-adjoint in this space. 

The use of nontrivial and, in general, manifestly Hamiltonian-dependent Hilbert-space- 
metric operators may be perceived as an important innovation of the model-building in 
quantum theory. Nevertheless, whenever accepted as a sound theoretical tool in physics, 
its mathematical consistency must always be carefully reexamined. In this sense the re- 
cent rigorous mathematical proof [12] of the nonexistence of the metric operator for the 
most popular and phenomenologically highly relevant H^^^\5) makes the study of this 
particular Hamiltonian far from being completed. For illustration purposes, fortunately, 
various other, mathematically simpler and fully consistent models are available and may 
be found in the current literature [5], | 



2.2 Bounded-operator models 

One of the main and well known [11] difficulties accompanying the use of the Hamiltoni- 
ans which are defined as differential operators lies in the necessity of a guarantee of the 
boundedness of the metric operator 6'^'^^ in Eq. ([T]). Many authors are trying to avoid 
this purely mathematical obstacle by the assumption of the boundedness of the operators 
of observables [3]. In this spirit we replaced, in paper [13], the above-mentioned ordinary 
differential Hamiltonian H^^^^^ by its truncated and shifted, A^— dimensional diagonal- 
matrix simplification 

""-(A^-l) ... 

-(iv-3) ■•. ; 
; ■•. ■•. 

... +(iV-l) 

In parallel, the methodical and pedagogical role of the Bender's and Boettcher's Hermiticity- 
violating differential-operator Hamiltonian H^^^\5) was transferred there to the multi- 
parametric anharmonic-oscillator-like real and tridiagonal non-Hermitian N hj N matrices 
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This allowed us to avoid the necessity of the complicated proof of the reality of the bound- 
state spectra of H^^^\5) (cf. [15]). For our models (|3]) we reduced the proof to a virtually 
elementary spectral-continuity or spectral- inertia argument, applicable at the not too large 
couplings Qj at least. Moreover, after an additional up-down symmetrization of the Hamil- 
tonian matrix, i.e., after the choice of Qn-i = Qi and gN-2 = 92, etc, we arrived at the 
final form of our benchmark toy-model matrices 
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(4) 



With such a choice we paralleled the phenomenologically relevant parity times time re- 
versal symmetry of the above-mentioned toy-model Hamiltonians H^^^^ (S) by a formally 
analogous PT— symmetry of our matrices (jl]). This merely required the specification of 
T as the transposition plus sign reversal. The parity-simulating indefinite square root of 
the unit matrix appeared then represented by the antidiagonal, by matrix V with 
elements Vj^N-j+i = 1, j = 1, 2, . . . , A^. 

In our subsequent paper [16] we turned attention to one of the methodically most 
welcome features of the "PT— symmetric and J = —parametric benchmark Hamil- 

tonian (jlj), viz., to the amazingly elementary geometric form of the boundary dV of the 
J— dimensional compact domain V of the real parameters Qj for which the spectrum of 
H^p^^ remains real. For our models (jl]) this boundary (or, in the language of physics, 
the horizon of the bound-state stability of the system [17]) has been shown to acquire, at 
any matrix dimension A^, the same generic geometric form of the surface of a smoothly 
deformed hypercube with protruded (hyper)edges and vertices (cf. also [T8| [T9] for more 
details). 

In paper I we finally restricted our attention to models (jl]) in the strongest-coupling 
extreme. In a sufficiently small vicinity of any one of the vertices we revealed that there 
exists a certain specific A— parametrization of the couplings 9j{X) such that 



(I) 

the two by two matrix H^prp^{\) appears useful as a benchmark model of an energy- 
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bifurcation QC scenario in which 

Eo = -VX, Ei = +VX, (5) 

i.e., in which the spectrum is real iff A > Aq = and in which the whole spectrum 
becomes completely degenerate iff A = while it finally gets purely imaginary iff 
A < 0; 

• the three by three matrix H^p^-^{\) may be used as a benchmark model of a new, 
energy-trifurcation quantum catastrophe in which 

Eq = -2v^, El = 0, E2 = +2v/A. (6) 

Again, the spectrum proves completely degenerate iff A = 0. Up to the exceptional, 
A— independent real level -E[Ar/2] = emerging at any odd A^, the rest of the spectrum 
is also purely real or imaginary iff A > or A < 0, respectively; 

• the four by four matrix H'^prp^{X) may then serve as a benchmark quantum model of 
the analogous energy- quadrifurcation catastrophe in which 

Eo = -3y/X, Ei = -y/\, E2 = y/\, ^3 = 371, (7) 

etc. Same conclusions were found to hold, for Eq. (jlj), at any integer N >2. 

In what follows we intend to develop further the emerging idea of an A^— furcating quantum- 
theory generalization of the classical bifurcation catastrophe. In this direction our present 
paper will clarify, first of all, some of the problems which emerge during the necessary 
reconstruction of the physical Hilbert-space metric G = 0(A). 

3 Quantum theory near the singularity 

3.1 The concept of the time of recovery t > 

One of the simplest forms of the above-mentioned A— parametrizations of the couplings in 
Hamiltonians @ is given by the formula which was proposed in paper [16], 

g^f^\\) = Vn{N-n){l-\)eV, n = l,2,...,N-l. (8) 
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The merit of this parametrization is that in the whole interval of A G (0, 1) (or for- 
mally even of A G (0, oo)) it guarantees the reality of the energy spectrum by keeping the 
whole J— plet of couplings safely inside the physical domain V. The second merit of the 
A— parametrization lies in the fact that the boundary value of A = strictly separates 
the dynamical quantum A > regime yielding the real, "observable" A^— plet of bound 
state energies from the completely unobservable A < half-axis. Geometrically this means 
that the value of A = localizes the protruded QC vertex of the domain V. In the algebraic 
terminology of original paper [T^ and of the older literature j20| [21] one encounters the so 
called extreme exceptional point (EEP) at A = 0. 

In paper we restricted our attention to the sufficiently small As. We emphasized 
that Eq. may be further reparametrized in terms of the time which would measure a 
recovery from the QC singularity, 

X Xn{t) =t + t^ + ... + t^-^ + Gnf^ n = l,2,...,Ar, J=[N/2]. (9) 

Formally, the time-parameter t G (0, oo) was redundant but its presence enabled us to 
treat Hamiltonian Hj^p^^ as manifestly time-dependent. 

3.2 Metric O near the singularity 

In the above picture the evolution is interpreted as the motion of the system away from 
QC, towards a more stable dynamical regime. Parametrization (|9]) proves useful, first of 
all, at the very short times t <^ 1 at which it effectively rescales and magnifies the interior 
of V in the EEP vicinity. Simultaneously, such an ad hoc change of scale did not lower the 
number of degrees of freedom - one could still work with as many as J coupling constants 

Gn > 0. 

With the trivial selection of the new coupling constants Gn = the presence of the 
time t allowed one for a one-parametric change of the dynamics. After the hypothetical 
start of the evolution at the QC singularity t = one reaches the manifestly Hermitian 
harmonic-oscillator limit at A(t) = 1. In this sense parameter A G (0, 1) may be perceived 
as an alternative, simplified measure of the QC-recovery time. 

At any fixed value of A^ and for any given Hamiltonian the physical contents of the 
theory varies with the inner product ([1]). The reconstruction of all of the eligible metrics 
B may be found summarized in our dedicated work [22]. It is based on a replacement of 
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matrix Hjp^^hj its transposition in Schrodinger equation, 



H, 



(N) ■ 
(PT) 



l^f^))=i^fMV^f^)), - = 0,1, 



(10) 



The complete solution of this equation opens the way towards the reconstruction of any 
metric from its spectral representation 

N 



fill 



n=l 



All of the parameters k„ > are freely variable. This, as we already mentioned, is an 
ambiguity which is usually being suppressed via a more detailed information about the 
system (cf., e.g., Ref. [3] for explanation). 



3.3 The ambiguity of the metric 

Let us pick up = 2 and study the possible consequences of the ambiguity of metric (ITTl) 
in more detail. Firstly, let us change the variables, t — )■ r = r(t) = \/t > 0, yielding 



H, 



(TV) ■ 
(PT) 



1 



(12) 



In terms of the pair of abbreviations u = a/1 — r and v = + r we may then follow 
paper I and calculate the maximal and minimal eigenvalues = r and = — r of (fT2|l 
as well as the related respective eigenvectors 



= [Vl-r, -Vl + r] = [u, -v] 



(13) 



|^_)) = [VTT^, -VT^] = [v, -u] . (14) 

It remains for us to insert vectors f|T3|) and f|T4|) in the spectral expansion of the metric. 
Once we fix an inessential overall factor and denote /t+ = sin a and k_ = cos a with 
0<a<7r/2we get the general N = 2 metric-operator matrix of paper 



1 + r cos 2a 



-VT 



— a/I — 1 — r cos 2a 
Its elementary form facilitates the direct determination of its eigenvalues. 



(15) 



e± = l±Vl- r^sm^2a. 



(16) 
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One easily verifies that the requirement of the necessary positivity of these eigenvalues is 
trivially satisfied at any square- root-time r = i/t such that < r < 1. 

We may conclude that the standard probabilistic interpretation of our time- dependent 
N = 2 THS QC quantum model is determined not only by its one-parametric Hamiltonian 
(fT2|) but also by the specification of the concrete value of variable a. Via Eq. (fT5|) this 
choice selects one of the eligible inner products ([1]). This makes the Hilbert space of states 
fully defined and unique. 

4 Minimal anisotropy constraint 

The concept of quantum catastrophes perceived as infinitesimally short processes is cer- 
tainly incomplete. Indeed, we saw that the full picture of evolution contains multiple 
Hamiltonian-independent parameters, sampled by a of Eq. f|T5|) . Thus, the theory based 
on the mere knowledge of the Hamiltonian would still admit many phenomenologically 
nonequivalent predictions. At non-infinitesimal times t the evolution would remain un- 
determined. Naturally, such a theory would be unsatisfactory. The ambiguity must be 
removed via additional, physics-based constraints [3]. 

One of the most popular options may be found advocated, e.g., in review It makes 
use of an additional postulate of the observability of a charge. In what follows, we intend 
to propose another, conceptually different form of the necessary constraint. We shall be 
able to make the metric unique via a new, QC-related postulate. 

4.1 A reversal of the arrow of time 

In paper I we were only able to construct the unique and closed-form metrics up to = 3 
or, in the implicit spectral-expansion form f ill I) , up to = 5. The puzzling arbitrariness 
of the choice of parameters k in Eq. (ITT]) was successfully eliminated by means of an ad hoc 
requirement of a maximal simplicity of the formulae for the matrix elements of the metrics. 
The main motivation of such a rather formal method of suppression of the unwanted 
ambiguity of the physical predictions was rather vague. The main argument supporting 
this strategy was that it perceivably facilitated the extrapolations of the formulae with 
respect to the dimension A^. 

The entirely artificial nature of such a half-intuitive and mo del- dependent choice of ks 
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appeared to be one of the strongest motivations of a renewal of our analysis of the problem 
of ambiguity. Our efforts climaxed in the successful resolution of the puzzle which will be 
described in what follows. Our new recipe will be based, first of all, on a removal of the 
restriction of attention just to a small vicinity of the QC instant t = 0. We shall extend 
the applicability of our toy models to the explicit description of the whole QC-generating 
process. 

The change of philosophy will force us to replace the not too suitable time-of-recovery 
t (or, if you wish, the time-like parameter X{t) of Eq. (j9])) by the new time variable r. It 
will run in opposite direction, i.e., from the r = instant at which our Hamiltonian was 
still entirely standard (i.e., in our toy models, diagonal) to the QC r = 1 instant at which, 
by construction, our Hamiltonians become non-acceptable - fully non-diagonalizable and 
merely Jordan-block representable |24] . 

4.2 Metrics as functions of r 

The choice of the concrete definition of the new "time of collapse" r = — A will enable 
us to treat the fall of our system into its singularity as a process which started at r = 0, 
i.e., long time before the catastrophe. A disadvantage of the replacement of t by r might be 
seen in the fact that the construction of the metric appears much easier at the small QC- 
distances A. Indeed, in the QC limit r — t- 1 the time-dependent vectors lipn)) are getting 
mutually parallel so that the operator ( |TT|) itself degenerates to a particularly elementary 
matrix of rank one. Fortunately, near the opposite extreme of r = the transition t — )■ r 
simplifies the Hamiltonians, 
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(17) 



In other words, one feels tempted to replace the explicit use of the spectral expansion 
fllip by some less subtle, brute-force construction techniques which make the use of the 
tridiagonal- matrix form of the Hamiltonians [25]. Such an idea proved truly productive. 
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Theorem 1 The metrics 6^^^(r) compatible with the respective Hamiltonians ( (771 ) may 
he sought in the generic form 

QW[r) = Y,{-ry-'M^'^\3) (18) 

containing some suitable sparse-matrix coefficients 

aii(l) ... 

ai2(l) ■■■ ■■ 
': ■■. ■■. 

... aiAr(l) 



(19) 



M(^)(2) 



aii(2) 
a2i(2) ai2(2) 







^22(2) ai3(2) 
... a2,7v-2(2) 
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a2i(3) ai2(3) 

«3l(3) ^22(3) 








ai,7v-i(2) 
a2,jv-i(2) 



ai3(3) 



(20) 



«32(3) 













a2,7v-3(3) «i,Af-2(3) 

a3,Ar_3(3) a2,7V-2(3) 

03,^-2(3) 



(21) 



etc. In this way, the set of all of the non-vanishing elements of matrix Ai^^\k) may be 
arranged into an auxiliary, k by {N — k + 1) — dimensional array 



a{k) 



aii{k) auik) ai3{k) ... ai^N-k+i{k) 
<y2i{k) a22{k) a23ik) ... a2,N^k+i{k) 



aki{k) ak2{k) aks^k) ... ak,N~k+i{k) 
at any k = 1,2, . . . , N , with aii{k) = M.[^\k), etc. 



(22) 
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Proof follows from the inspection of the set of the A^^ Dieudonnes linear algebraic com- 
patibility relations written in the matrix form 

H^Q = QH (23) 

(not all of which are independent - cf. Ref. [22] for details). □ 



4.3 The physics of the onset of degeneracy 

From the point of view of potential applications of formula f fTSj) it is important that at the 
beginning of the fall into instability (i.e., at the far-from-QC instant r = 0) the Hamil- 
tonians f[T7l) will all coincide with the respective truncated and diagonal (i.e., Hermitian) 
harmonic-oscillator-like matrices. In this picture the spectrum of energies Ei''\T) will 
remain real but shrinking with the growth of the innovated time r. In the limit r — 1, 
i.e., at the very end of the fall of the system into QC singularity, the spectrum becomes 
completely degenerate, Eif^\l) = 0, n = 0, 1, , . . . , — 1. 




0.2 0.4 0.6 0.8 1 

t 



Figure 1: Eigenvalues ( [T6l) of the N = 2 metric ( fT5i) for the recovery time t = r"^ & (0, 1) 
and for the parameters (a) sin^ 2a = 3/5, (b) sin^ 2a = 4/5 and (c) sin^ 2a = 1. The 
desirable disappearance of the anisotropy at the full-recovery time t = 1 only occurs for 
the minimal-anisotropy choice of sin^ 2a = 1. 

At the latter instant the Hamiltonian H^^\l) ceases to be diagonalizable and loses 
its standard physical tractability and interpretation. Vice versa, from the point of view 
of physics the description of the evolution as generated by H^^\t) will change at r = 1, 
requiring an introduction of some new degrees of freedom beyond this instant, i.e., at 
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r > 1. The study of such a switch to a new form of Hamiltonian at later times hes already 
beyond the scope of the present paper. Examples of such a scenario may be found, e.g., in 
magnetohydrodynamics [IB] . 

Temporarily, let us now return, in the context of interpretations, to the times of recovery 
t or A(t). In these variables, the motion beyond the end-of-the-interval A(t) = 1 appears 
much less exotic. Typically, the energies would not feel the change at all (cf., e.g., Eqs. ([5]) 
- dZD). Still, the values of A = A(°"*^^) > 1 remain mathematically exceptional because in 
the "outer" interval of parameters our Hamiltonian matrices change their form thoroughly: 
they become complex and Hermitian. 

We are now coming very close to the key message of our present paper, referring to 
the fact that in the "outer" interval of parameters our Hamiltonians become manifestly 
Hermitian in usual sense. For this reason it seems most natural to require that the metric 
remains trivial in such an interval, Q(°"*^^) = J. In other words, we suggest to match the 
separate non-Hermitian and Hermitian dynamical regimes strictly at the natural matching 
point with X(t) = 1. 

One of the consequences of such an innovative model-building decision may be seen in 
the subsequent most natural reinterpretation of the point of matching A(t) = 1 as the time 
of onset of the whole QC scenario using A(t) G (0, 1). The optimal QC-related metrics 
should be then required continuous at the instant X{t) = 1, i.e., at the very onset of the 
whole process of the QC degeneracy, 

lim e^^Hr) = 1. (24) 

T— >0 

In a climax of this paragraph we must add and emphasize, incidentally, that the new 
restriction proves strictly satisfied by the family of the formally admissible metrics in 
their unique extrapolation-friendly form as introduced in paper I. 

The latter observation may be most immediately illustrated via the N = 2 model 
again. Recalling the general set of the metrics f|T5|) which appear numbered by the sin- 
gle kinematical real variable a, we may conclude that each deviation of the latter free 
parameter from its extrapolation friendly value of paper I and from its anisotropy min- 
imizing value which follows from formula (fT6l) also, simultaneously, violates the present 
onset-continuity constraint (12^ . In this sense, all of the three forms of the constraint 
may be perceived equivalent. For general dimension A^, their most immediate test will be 
provided by Eq. of course. 
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Figure [T] offers the simplest N = 2 illustration of the fine-tuning nature of the coinci- 
dence. We see there that the influence of a changes from very weak in the QC vicinity (i.e., 
at t ^ 1) to very strong near the Hermitian LHO dynamical regime (where X{t) = 1). 
In such a graphical representation we also see that the extrapolation-friendly choice of 
a = n/A oi paper I is truly exceptional. Solely for this choice of the free parameter the 
difference between the eigenvalues of the metric will strictly vanish at X{t) = 1. We may 
just repeat that for the value of a = 7r/4 the r = instant really carries the meaning of a 
Hermitian onset of the fall into QC singularity. 

Once we preserve the latter exceptional choice of the parameter a = 7r/4 at all times 
t G (0, 1), the difference between the two eigenvalues of the N = 2 metric (measuring a 
Hilbert-space anisotropy) remains always minimized (cf. both Fig. [T] and Eq. flTBl) ). In 
this formulation, the selection of a "minimal Hilbert-space anisotropy" principle certainly 
becomes much less model-dependent. Moreover, the uniqueness of the N = 2 choice of 
a = tt/4 may be, mutatis mutandis, extended to all of the higher matrix dimensions A^. 
After all, one might most easily proceed from the point of the onset r = 0, i.e., from the 
very end A = 1 of the recovery of the system from the QC dynamical regime, by connecting 
our specific removal of the mathematical ambiguity of the metric directly with the well 
founded physical trivalization (12^ of the metric in the Hermitian-Hamiltonian limit. 

5 Hilbert spaces with minimal anisotropy 

In paper [22] we thoroughly discussed the general recipe by which suitable metrics may 
always be assigned to a given Hamiltonian H via Eq. fITT]) . We reemphasized there that 
such a simple-minded construction is always ambiguous [3] but that for the finite space- 
dimensions N < oo the key merit of the recipe lies in the sufficiency of the solution of the 
auxiliary conjugate Schrodinger equation (ITU]) . 

In paper I we applied this philosophy to Hamiltonians (|T71) . Up to = 5, we managed 
to evaluate the necessary input ketket-eigenvectors | tpi^"^ )) . We also revealed a very regular 
sparse-matrix pattern in these formulae [cf. Eq. Nr. (10) in paper I]. Finally, we empha- 
sized the universality of the recipe flTT]) yielding all the admissible metrics. At the same 
time, our interest remained restricted to the mere vicinity of the QC dynamical regime. 

In our present paper we shifted the emphasis from the universality of spectral for- 
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mula f|TT]) to the alternative possibility of its explicit resummation. Let us now test this 
strategy at = 2 once more. 



5.1 N = 2, revisited 



The explicit and universal construction of the metric Q^^\t) as presented in paper I (i.e., 
the one performed via the auxiliary Schrodinger Eq. ( fTOj) ) is not too easy even at iV = 2, 
i.e., for our first nontrivial QC-related Hamiltonian matrix 



T 



-1 r 
-r 1 



(25) 



This may make the efficiency of this construction comparable with the direct solution of 
Eq. fl23l) . For a verification, let us now return to the latter methodical possibility and let 
us try to rederive, say, the complete family of the N = 2 metrics 0*-^^ as already known 
from Ref. [23] and also from subsection 13.31 above. 

In the real-matrix ansatz with the subscripted vector k containing two arbitrary positive 
components, 

a b 



0gW 



b d 



(26) 



we shall, first of all, fix an overall multiplication constant by setting the determinant equal 
to one. This enables us to put b = sinh u and choose e = ±1 in a = e cosh u exp g and 
d = e coshz/ exp{—g) (both the new parameters u and g are assumed real). As long as 
the metric must be positive we may only use e = 1. Finally we check that the matrix 
constraint ( 12^ degenerates to the single, time-reparametrization item 

tanh u 



T 



(27) 



cosh g 

Our conclusion is that for any given r G (0, 1) we may choose any real g G (0, gmax) (note 
that this is the parameter which makes the main diagonal of the metric asymmetric). This 
choice enables us to evaluate u = z/(r, g) from the latter equation (this implies that at 
a fixed time, the value of gmax must be such that cosh gmax = V'^)- Summarizing, we 
may set aii(l) = coshu exp g, auil) = cosh z/ exp (—f?) and aii(2) = sinhz/ in Eq. f lTSj) at 
N = 2. The resulting eigenvalues of the metric are both, by construction, positive. 



9± = cosh 1/ cosh ± y cosh z/ cosh g — 1 



(28) 
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At the very start of the fall of the system into the catastrophe, i.e., at r = one has 
Pmax(O) = oo so that there is no upper bound imposed upon g{0). Still, as long as one 
might like to have the trivial, isotropic initial value of O^^^(O) ~ / (implying the special 
choice of u{0) = and g{0) = 0) the resulting metric becomes, up to the above-mentioned 
irrelevant overall multiplication factor, unique at r = 0. 

During the subsequent growth of r < 1 the requirement of the minimization of the 
anisotropy leads to the rule ^{t) = (cf. Eq. (1281) ) so that the remaining variable u < 
may now be interpreted as another (viz., rescaled and, incidentally, inverted, cf. Eq. fl27j) ) 
version of the time of the QC degeneracy. 

Once we return to the standard variables we get our unique and minimally anisotropic 
metric in the virtually trivial form as given in paper I, 

1 -T 
- T 1 

From this formula we may deduce the special, minimally anisotropic version of eigenvalues 
in the form compatible with their "more anisotropic" generalization f|T6l) . 



Q(2) 



I-tJ. 



5.2 N = 3 

Whenever one tries to move to the higher matrix dimensions one encounters the technical 
problem of the increasing multitude of various parameters. Thus, in the first nontrivial 
case with = 3 let us first follow the N = 2 guidance (cf. the ultimate choice of g = 
in the preceding paragraph 15. ip and let us omit the discussion of the metrics with an 
asymmetric form of their main diagonal. 

Once we also keep ignoring the other, irrelevant though still existing overall factor, 
we are, after some straightforward manipulations using Eq. (1251) . left with the last free 
parameter g in the metric 

1 -V2gT 

-V2gr 2g-l + gr^ 

gr^ -V2gT 
Among its three readily obtainable eigenvalues 

ei = gT'+g-^Ag^T^ + g^-2g + l, 



6(3) (r) 



gr^ 

'V2gr 
1 



(29) 



9. 



l-gr' 



9, 



gr'^+g+^Ag^T'^ + g'^ -2g + l 

(30) 
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the middle one (i. e., the inverted parabola in r) remains positive for the parameters 
g < while the change of the sign of the remaining pair may take place at the curves 
g = and g = l/{2 — r^) in the g — r plane. Figures [2] - S] may be recalled as showing 
why the correct and unique choice of the parameter is g = 1, yielding again the unique 
metric of paper I, 

72 
V2 

with the correct and expected r— dependence of the eigenvalues as given by Eq. fl30|l . 



(31) 
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Figure 2: Eigenvalues ( l30l) of the = 3 metric (l29i) in an extended range of time of fall r 
and for the "correct" parameter g = 1. 
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Figure 3: Eigenvalues ( 130|) of the = 3 metric (1291) in an extended range of time of fall 
r and for a "too large" parameter g = 1.2. Notice that the metric is not isotropic even at 
the onset time r = 0. 
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Figure 4: Eigenvalues f l30|) of the = 3 metric fl29!) in an extended range of time of fall r 
and for a "too small" parameter g = 0.8. Notice that in this regime the metric ceases to 
be acceptable (i.e., positive) long before the QC time r = 1 is reached. 

5.3 iV = 4 
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Figure 5: Eigenvalues ( l33l) of the "correct" = 4 metric ( l32i) in an extended range of 
time of fall r. 



In the next step of our constructive considerations we are getting beyond the formulae 
derived in paper I. Fortunately, already the = 3 formula (IM]) itself offers a clear hint of 
an extrapolation so that it just proves sufficient to verify that the next metric 



e(4) 



1 



V3i 



VSt 1 + 2 -2 r - ^3 



r 



V^7 



-V3r 

-r- V3t^ -V3t 1 
obeys all the necessary and sufficient requirements including, first of all, the Dieudonne 



-2 r - 1 + 2 r2 
-2 -v^r 



(32) 
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Eq. fl23l) . Once we also deduce the r— dependent eigenvalues of this candidate for the 
metric at = 4 we may immediately see that they behave as they should, 

{^1, . . . , ^4} = {1 - 3 r + 3 - r^ 1 - r - + r^, 1 + 3 r + 3 + r^, 1 + r - - r^} 

(33) 

(cf. also Fig. [5] for a graphical illustration of their correct QC behaviour at r = 1). 

6 Extrapolations in N 

6.1 Metrics between = 5 and N = 7 

We may now combine the results of preceding section with the contents of Theorem 1. 
Using an elementary insertion in Eq. ( p3|) we may easily prove that in the expansions 
( fTSl) of the metrics with minimal anisotropy the diagonal-matrix coefficients ( !T9l) may be 
defined, at all A^, by the elementary formula 

ai„(l) = l, n = l,2,...,N. 

Similarly, the closed formula is also available for the antidiagonal coefficients in A^^^)(A^), 

a„i(iV) = l, n = l,2,...,N. 

Next, the bidiagonal matrix coefficients (l20l) may be defined, at all A^, by the slightly less 
elementary general formula 

ain(2) = a2n{2) = ^/n{N - n) , n = 1, 2, . . . , iV - 1 . 

Due to the easily verified symmetry, the analogous formula exists for the coefficients in 
yuW(iV-l), 

«„i(Ar-l) = «„2(iV-l) = V^iN - n) , n = l,2,...,Ar-l. 

Up to now, unfortunately, we did not succeed in an extension of these observations to the 
tridiagonal sparse matrix coefficients ( 12T|) . etc. Nevertheless, we believe that the task is 
not impossible. This belief seems supported by Theorem 1, i.e., by the reducibility of the 
N by N matrices Ai^^\k) with A; = 3,4, ... to the respective auxiliary k hy N — k + 1 
arrays containing the non- vanishing matrix elements ajm{k) of Ai^^\k). 
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The first missing set of coefficients occurs at = 5. Its values 
aii(3) = ai3(3) = a3i(3) = 033(8) = Vg , 

"12(3) = 021(3) = 023(3) = 032(3) = 3 , 022(3) = 4 . 

may still be found evaluated, albeit in different context, in paper I. Naturally, this definition 
should be better rewritten in the much more compact form of the array 



0(3) 



^6 3 ^6 
3 4 3 
76 3 V6 



N = 5. 



(34) 



It makes sense to complemented this result by the next, = 6 formula 

Oii(3) = 014(3) = 031(3) = 034(3) = v^, 021(3) = 024(3) = 4, 

012(3) = 013(3) = 032(3) = 033(3) = 3 V2 , 022(3) = 023(3) = 6 

which we derived using the brute force construction based on Eq. (j23|) . It again deserves 
the compact presentation as the array 



0(3) 



To 3V2 3V2 Vm 

4 6 6 4 

To 3V2 3V2 VTo 



N = 6 



(35) 



The closed form of the latter result indicates that there might exist a not too complicated 
extrapolation recipe, with the help of which we would be able to determine the unique, 
minimally anisotropic metric at any dimension A^. This belief seems further supported 
by the regularity and apparent extrapolation-friendliness of the next two sparse-matrix 
"missing" coefficients 



A^(^)(3) 









V15 

















5 
















/T5 





8 





6 











V30 





9 





V30 











6 





8 





















5 

















^Ts 
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and 


















2 








2v^ 


6^3 


2v^ 




2^5 


6^3 


6^3 






2v^ 


6^3 









2v^ 


2 
















Again, they were obtained, with the assistance of the computerized symbohc manipula- 
tions, by the brute-force solution of the set of 49 linear algebraic Eqs. fl23l) . 



6.2 Eigenvalues of the metrics 



Table 1: Pascal triangle for coefficients C\J in Eq. ( I36ll 
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10 




5 




1 


7 
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15 




20 




15 
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1 


8 


1 


7 




21 




35 




35 




21 




7 1 



In the above-described constructions of the N hj N matrices of metric O^^^^ we did not 
manage to find, unfortunately, any obvious general extrapolation tendency or pattern. For 
this reason, we turned our attention from the matrices to the perceivably simpler-to-display 
A^— plets of their eigenvalues 6^\t). 
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Table 2: Pascal-like triangle for coefficients in Eq. f l36|) 
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5 
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5 




-5 




-9 




-5 -1 



Proposition 2 The time- dependent eigenvalues of Q^^^ may be written in the fi 



orm 



N 



0L''^ = 9ir\r) = J2ciVr'-' (36) 

k=l 

where, up to N = 8 , the numerically evaluated values of the coefficients C^^^ may be found 
summarized in the Pascal-like tablesU\-^ 

In the final step of our considerations let us now show that the closed-form eigenvalues 
of the metric can be obtained at any matrix dimension A^. 

Theorem 3 The time- dependent eigenvalues of Q^^^ are given by formula 

N 

0f\r) = Y.Ci'2r-\ k = l,2,...,N 

m=l 

where Cl^ = \ \ , =1 ~ I I '^'^^^ general, 

n-1 I \ n-1 j \ n-2 

Proof is straightforward and proceeds by mathematical induction. 

□ 





II 


( N-k\ 






\ n — p j 
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Table 3: Pascal-like triangle for coefficients C^^ in Eq. f l36|) 
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Table 4: Pascal-like triangle for coefficients C^^ in Eq. ( !36ll 
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7 Discussion 

We may summarize that our present, toy- model-based generic scenario of the evolution 
of a quantum system S sets the scene between its innocent-looking and safely Hermitian 
equidistant-energy-level onset prepared at the initial time-of-fall r = and a full realization 
of the complete, A^— tuple degeneracy of the energy levels at the final QC time r = 1. 

According to the general principles of quantum theory in its THS form the states of the 
system S are assumed described by a wave function tp which evolves in time in a way which 
has been thoroughly described in Ref. [6]. In our present paper we omitted the repetition 
of the related discussion and never mentioned this purely technical and routine aspect 
of the model. We rather restricted our attention to the details of the interplay between 
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the time-dependence of the pre-selected benchmark Hamiltonians H^^^ and of the exactly 
constructed and minimally anisotropic unique metric operators Q^^\ 

The set of our conclusions offered a compact and consistent picture of the process in 
which the not quite expected exact solvability of our toy model of paper I enabled us to 
extend its applicability to all times r G (0, 1) measuring the whole QC-preceding history 
which started, at r = 0, from an entirely standard, Hermititian- matrix equidistant-level 
initial Hamiltonian. Via our A^— numbered family of benchmark models given and defined 
by the matrix pairs 

we were able to explain the quantum-evolution fall in the benchmark, A^— tuple level- 
degeneracy catastrophe as a consequence of the steady growth of the anisotropy of the 
underlying physical, manifestly time-dependent Hilbert space 'H^'^^ 

We found it natural to characterize the extent of the latter Hilbert-space anisotropy 
by the non-coincidence and certain "spread" of the set of the eigenvalues of the pre- 
selected physical metric Q^^\ On this background we decided to make the metric unique 
by means of the requirement of the minimization of this spread at any time r and, more 
explicitly, by the ultimate and decisive requirement of the zero limit of the special spread 
measure p = ma,x{6li^^ — 6^^) at the onset r = of the whole hypothetical evolution 
process. On this background we finally arrived at the family of benchmark QC models in 
which our exceptional and unique metric 6 remained nontrivial during all times r G (0, 1). 
At the beginning of the QC history with r = 0, this metric was assumed to be matched to 
the trivial and most common isotropic metric ©('D*'''^^) = /. In such a manner we arrived 
at the benchmark quantum representation of the catastrophe in which the fall into the 
A^— tuple degeneracy as running in the interval of time r G (0, 1) may be perceived as 
initiated and smoothly "switched on", at r = 0, to its arbitrary hypothetical unitary- 
evolution r < prehistory. 
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